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Abstract 

We consider a reaction-diffusion-advection equation arising from a bi- 
ological model of migrating species. The qualitative properties of the 
globally attracting solution are studied and in some cases the limiting 
profile is determined. In particular, a conjecture of Cantrell, Cosner 
and Lou on concentration phenomena is resolved under mild condi- 
tions. Applications to a related parabolic competition system is also 
discussed. 
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1 Introduction 

In mathematical ecology, reaction-diffusion equations are often used to deter- 
mine the factors behind the survival and extinction of animal populations. 
(See for examples [H El El II])- One well-known example is the following 
logistic reaction-diffusion model for population dynamics (See [5]): 
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Ut = dAu + u[m{x) — u] in Q x (0, oo), 
1^ = on Onx (0,oo), 



(1) 



where u{x, t) represents the population density, A = Ylii=i ^ Laplace 

operator in R^, d > is the dispersal rate, m{x) accounts for the local growth 
rate, VL is the habitat of the population and is assumed to be a bounded re- 
gion of with smooth boundary dVt, and v is the outward unit normal 
vector on dVt. The Neumann boundary condition, which coincides with the 
no-flux boundary condition, is imposed on dVt. 

If the environment is spatially heterogeneous, i.e. m{x) is non-constant, then 
it seems reasonable to assume that the population has a tendency to move 
up the gradient of m{x) in addition to random dispersal. In this direction, 
Belgacem and Cosner [6] proposed the following reaction-diffusion-advection 
equation: 



where the parameter a > measures the rate at which the population moves 
up the gradient of m{x). Again, the corresponding no-fiux boundary condi- 
tion, is imposed. For discussions on the modeling aspects, we refer to [6l [7] 
and the references therein. 

The dynamics of ([2]) seems simple. In fact, it was established in [21 [S] that 

if we assume that 

(HI) m(x) G C^{VL), and is positive somewhere, 

then for any d > ^ has a unique positive steady-state u for all large a. 
Moreover, u is globally asymptotically stable among all nonnegative, nonzero 
solutions. In other words, the steady-state u of ([2]) determines the long-time 
behavior of all solutions of ([2]). We shall always assume (HI) throughout 
this paper. 

From both mathematical and biological points of view, it seems important 
to understand the qualitative properties of u. In particular, it would be 
interesting to describe the shape of u. There has been considerable effort in 
this direction. Recently, it was proved in [S] that if the set of critical points 
of m{x) has Lebesgue measure zero, then 



Mt = V ■ {dWu — auWm) + u{m — u) in Q x (0, oo), 
d^ - au^ = on dnx (0, oo) 



(2) 
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That is, the total population size tends to despite the fact that the species 
is tracking the resources more accurately. To understand the mechanism 
behind such phenomenon, again a better description of the shape of u is 
desired. To this end, the following results were proved. 

Theorem 1.1 (Cantrell-Cosner-Lou). Suppose m{x) > in Q. Let u he the 

unique positive steady-state of ([2D. 

(i) If a > d/ min^-m, then u{x) > max^ym ■ g"("i(^)-ma,Xfj-m)^^ every 
X & Q. In particular, max^u > max^ym. 

(a) SupposeQ = (—1,1), andm{x) has finitely many critical points {xi}f^i, 
then M — i- uniformly in compact subsets ofQ \ {xi}^^i as a oo. 

Based on these results, the following conjecture was proposed in [9] and 
Section 3.2 in |TD]. 

Conjecture 1.2. u concentrates precisely on the set of (positive) local max- 
imum points of m{x) as a ^ oo. 

Remark 1.3. We have modified the concentration set to be the set of positive 
local maximum points instead of local maximum points stated in fS] , since we 
are considering a more general situation where m{x) can change sign on the 
set of its local maximum points. 

In this paper we shall establish Conj ect ure 11.21 under mild conditions on m{x). 

Let OJl be the set of all positive strict local maximum points of m{x) (i.e. 
those lying in {x & Q : m{x) > 0}). 

Theorem 1.4. Assume that u is the unique positive steady-state o/ ([2]). // 
xq e Wl, then for any ball B centered at xq, 

liminf sup u > m(xo). (3) 

In other words, u concentrates at each point of 971. The proof of Theorem 11.41 
is based on the observation that u solves a corresponding eigenvalue problem 
and is given in Section 2. 

To prove that u concentrates precisely on 9Jl, we impose the following as- 
sumptions on m{x). 

(H2) ^ < on an. 
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(H3) m{x) has finitely many local maximum points in Q, all being strict 
local maxima located in the interior of Q. 

(H4) Am(xo) > if Xq € is a local minimum or a saddle point of m{x). 

Theorem 1.5. Assume m{x) satisfies (H2), (H3) and (H4), then for any 
compact subset K of Q \ DJl, there exists 7 = '~f{K) > 0, such that 

< u{x) < e-^°, for all x e K. 

In particular, u — uniformly and exponentially in K, as a —t- 00. 

Theorems 11.41 and 11.51 together guarantee that u concentrates precisely on 
971, the set of positive local maximum points of m{x), thereby Conjecture 11.21 
is established. Theorem 11.51 is proved in Section [2] by the construction of an 
upper solution closely related to the shape of m{x). 

The question of determining the profile of u is, however, far more challenging. 
We only have the following result by a very interesting method introduced in 
[H] for the special case when m{x) is constant on the set of local maximum 
points of m{x). 

Theorem 1.6. If m{x) satisfies (H2), (H3) and (H4) and moreover, 

det -D^m(xo) 7^ for all Xq E 9Tt, 
with m{xo) = mi > for all local maximum points Xq G Q, then 

lim II u{x) - 2^/2rriie°[^'"(^)-'^i]/^ ||L-(m= 0. (4) 

a— !>oo 

Remark 1.7. The factor though mysterious at first glance, is ac- 

tually the consequence of the profile of u at each of its "weights", which is 
like a Gaussian distribution q°'[{^-^o)'^ D'^rn{xo){x-xo)]/2d ^ integral 

constraint j^i^^^^ — umdx = 0{e~'^°') for each xq G VJt. 

As in [9], [12], our resolution of Conjecture 11.21 has implications for the fol- 
lowing competition system. 

r Ut = V- {diVU - aUVm) + U{m -U -V) m ^] x (0, 00), 

I Vt = d2AV + V{m-U -V) m^]x(0,cx)), (5) 

[ di% - aU^ = 1^ = on X (0, 00). 

This system was introduced to model the competition of two species whose 
population densities are denoted by ?7(x, t) and V{x^ t) respectively. The two 
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species have identical local growth rate m(a;) and competition abilities, but 
different dispersal strategies: the species with density V disperses randomly, 
whereas the other species U disperses, in addition to random diffusion, by 
a directed movement towards more favorable locations, i.e. where m{x) is 
large. The goal of this model is to understand how different dispersal strate- 
gies affect the outcome of the competition in a heterogeneous environment. 

When a = 0, it is well-known [13] that if di > d2, then ([5]) has no coexistence 
steady- states, and solution {Ua.Va) o/© always converges to (O,^^^ o^st ^ 
oo, where is the unique positive solution to 

( d2A9 + 9{m-9) =0 in Q, 

1 ^ =0 on on. 

However, for any di,d2 > 0, the existence of the positive steady-states 
UayVa > of was established in [9l fTI] for all large values of a. Moreover, 
they proved that at least one of the co-existence steady-sates is stable! Some 
qualitative properties of these co-existence steady-states were also obtained 
under extra hypotheses on m{x). 

Theorem 1.8 (Chen-Lou). Suppose that J^m{x)dx > and all critical 
points of m are non- degenerate (detZ)^m(xo) 7^ 0). Then for any positive 
steady-state {Ua.Va) 0/®, 

liminf maxf/^ > max[m — 9dJ\ > 0, 
where 9d2 is the unique positive solution to . 

Assume further that m{x) satisfies (H2) and that m{x) has exactly one crit- 
ical point Xq which is a non- degenerate local maximum in the interior ofQ, 
then for any positive steady-state {Ua, Va) of 

V/3 e (0, 1) : lim II Va - 9d2 \\c^+P(n)= 0, and 
hm II f/^(a;)e-[--o"(-)]/'^i _ 2^/2[m(xo) - 9d,{x^)] \\L^in)= 0. 

a— s>oo 

Note that the condition J^m{x)dx > is there to ensure the existence of 
9d2- (See [9j.) It is interesting that our methods for ([2]) can be applied to 
study the coexistence steady-states. 

Theorem 1.9. Assume J^m{x)dx > 0. 
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(i) Assume that (H3) holds. Given any positive steady- state {Ua,Va) of 
O, if Xq G Tl, then for any ball B centered at Xq, 

liminf sup Ua > rn{xo) — Oa^i^xo). (7) 

a—^oo 

If in addition, (H2) and (H4) hold, then, for each compact subset K 
of fl\ VJt, there exists a constant 7 = 'y{K) > such that whenever 
{Ua, Va) is a positive steady-state of @, 

Ua{x) < e^^" for every x E K. 

(11) If (H2), (H3) and (H4) hold, deW'^m{xo) ^ for all xq G Tl, and 

m{xo) = mi > for all local maximum points xq G Q, then 

lim II K - e,, ||ci+^(f^)= V/3 G (0, 1), (8) 

lim II Ua{x) - 2^/\rm - (xo))e"["^(^)-™^l/'^i ||l-(oo= 0, (9) 

where Oi is any open neighborhood of xq such that xq ^ Oj for any 
other Xq G OJt. 

Remark 1.10. (i) ([7]) is useful only when m(xo) > 61(^2 (xq). And this is 
true on DJl if d2 > is sufficiently small and Am(xo) > 0. (The proof 
of this fact is included in Appendix A.) 

(a) The choice of 'y in Part (i) of Theorem \1.9\ is independent of choice of 
positive steady-state {Ua,Va)- 

(Hi) By maximum principle, mi — (2^0) > m Q for any d > 0. 

The rest of the paper are organized as follows. In Section [2] we provide the 
proofs for Theorems 11.41 II. 5[ and II. 61 Section [3] will be devoted to proving 
Theorem 11.91 Finally, some concluding remarks will be included in Section 

SI 



2 Proofs of Theorems II. 4L II. 5L and 11.6 



To simplify the presentation, we set = 1 in the proofs. This assumption can 
be removed with minor corrections. We first obtain the following equation 
for u: 



V ■ (Vm — auVm) + u{m 
^ - au^ = 

au av 



U 







in VL, 
on dVL. 



(10) 
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Proof of Theorem 1.^. Let u be the unique solution to f ITU]) . and Xq be a 
strict local maximum of m{x). Then u is the principal eigenfunction of the 
following eigenvalue problem with principal eigenvalue 0: 

V ■ (V0 - acjNm) + (m - u)0 + = m 



1^ - = on dVL. 



11) 



Now by the transformation (p = e"'"^ip, f jTTj) is equivalent to 

V ■ (e"™VV') + (m - M)?/)e""' + Ae"™'^ = m fi, 
1^ = on dn. 

av 

with principal eigenvalue equal to 0. The variational characterization of the 
principal eigenvalue of (|T2|) implies 



(12) 



= A = inf 



Given any small ball B = Brg{xo) centered at xq, since m{x) attains a strict 
maximum at Xq, maxoBrgixo) ''^^ < '"^(^^o)- For any e such that < e < 
— maxg^ (^.g) T7i, define 

e 2e 
Ml :=m(xo) " 3 > "^(a^o) - y := ^2, 

f/i :={x G 5r(,(a;o) : m(a;) > m{xo) - ^} 

t/2 :={x G 5,,o(a;o) : "^(x) > m(xo) - y} 

[/s :={x G (a^o) ^ "^(a;) > m{xo) — e}. 

Note that we have Ui CC U2 CC t/3 CC -B^ola^o)- Now take a smooth test 
function ip such that, 

= ^^^n\t/3 o^*Wsi WI<CW 



Then, 
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4 e"^i 4 6^-7^2 

< C"(e)e"(^^2-^^^) + niax(M - m) 

Us 

< C'{e)e~~ + rnaxM — m(xo) + e. 

c/3 
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For OL sufficiently large, the first term in the last line will become less than 
e, hence ([3]) follows. 

□ 

Next, we turn to the proof of Theorem ll.Si We first give the following 
definition of an upper solution. Denote from now on 

L0 = V ■ (V0 - Q!0Vm) + (m - 0)0. 

Definition 2.1. u is said to be an upper solution of f 1 1 01) if (i) (Hi) below 
hold: 

(i) There exists an open cover {Ui\ of Q, i.e. Q = [jUi where Ui's are 
relatively open in Q, and, 0j G C'^{Ui), Lipi < 0, such that 

u = min{0j} is continuous in Q. 

i 

(a) Denote Qi = {x & fl :u = 0j}. dQi is piecewise , and 

Qi CC Ui for all i. (13) 

(Hi) H — > for any x G dVl, whenever the normal derivative || is 

defined. 

The definition of lower solution can be obtained by reversing all the inequal- 
ities above and replacing min by max. 

The following is the key to obtaining an upper bound of u. 

Lemma 2.2. Fix a sufficiently large so that the unique positive solution u 
of (|T0|) exists. Ifu > is an upper solution of ( ITOj) in the sense of Definition 
l2l\ then u>u. 

To prove Lemma \2.2\ we first relate the above definition of upper solution to 
that of a weak upper solution from |14j . 



Definition 2.3. u G Vr^'^(f2) is said to be a weak upper solution of (ITU]) if 
it satisfies 

j ~ auVm) ■ Vip + u{m - u)il)\ < 0, for any ip G ly^'^(fi), ^/^ > 

1 1^ - au^ > on 

The definition of weak lower solution can be obtained by reversing the in- 
equalities appropriately. Note that by (H2), — > on dVL. 
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The following lemma can be proved via integration by parts. 

Lemma 2.4. Suppose u is an upper solution of O0|) in the sense of definition 
\2.1\ then it is a weak upper solution of f fTOj) . 

Remark 2.5. Lemma \2^ is true even if we drop the regularity of dfli in 
Definition \2.1\ provided we use the arguments in Lemma 4-10 of f75|. This 
observation will not he used in this paper. 

We recall the following well-known theorem on upper and lower solutions. 

Theorem 2.6 (Sattinger). Ifu and u are weak upper and lower solutions of 
fl 1 op respectively, and u > u, then there exists a classical solution u of OOjl 
such that u<u<u. Moreover, u is stable from above. 

We can now prove Lemma 12.21 by making use of the dynamics of ([2]) . 

Proof of Lemma \2.S[ Since u and are weak upper and lower solutions of 
fllOp respectively. By Theorem 12. 6[ there exists a solution u' which is stable 
from above such that < u' < u. Since is unstable in f lTOj) (by the global 
stability of u), u' ^ 0. Hence, u' = u (by the uniqueness of u). Therefore, 
we have u <u. □ 

To prove Theorem 11.51 it remains to construct an appropriate upper solution 
of f lTOj) according to Definition 12.11 To avoid complicated notations and to 
illustrate the ideas more clearly, we shall only prove in detail the cases: 

(a) When m{x) = mi > on 97t and m > at each of its critical points, 

(b) When m{x) = mi > on SOT and m < at some of its critical points, 

(c) When m{x) has two distinct values < mi < m2 on VJt and m < at 

some of its critical points. 

We remark that the same technique can be applied to prove the general case 
when m{x) has any (finite) number of distinct values on 071. The precise 
statement of the lemma that leads to Theorem 11.51 and some comments on 
its proof are included in the Appendix B. 

Proof of Theorem \1.5[ Case (a): When m{x) = mi > on OJt and m > at 
each of its critical points. 
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Lemma 2.7. Suppose that mix) satisfies (H2), (H3) and (H4). Assume 
m{x) = nil on SOT and m > at each of its critical points. Then for any 
a < 1, sufficiently close to 1, and for any < e < 1, there exists ao{e, c) > 
such that 

^€a{Tn{x)—CTni) 



ui = e 



is an upper solution of i \10^ in the sense of definition \2. 1\ for all a > ao. 
Proof. 

Lui = Aui — aVm ■ Vui + (m — mi — aAm)ui 

= ui{{e^ - e)«2|v^|2 ^ _ i)c,Am + m - e^"("^-'="^i)} 
= IIi{(e - l)a[ea|Vm|2 + Am] + m - e^°("-^"i)}. 

It suffices now to prove that the sum in the large parenthesis is negative. 

In {x G : m(x) < c^rrii}, by (H4), there exists ki > such that 

ea;|Vmp + Am > ki for all a large. 

While m — e'^"(''"~'=™'i) is bounded from above by |m|oo, therefore Lui < for 
all a sufficiently large. 

In {x e n : m{x) > y/crrii}, e'"("'-'="'i) > e'^"(v^""i-^"i) = e^^^ for gome 
^2 > 0. Whereas (e — l)a[ea|Vmp + Am] + m grows at most in the order 
a^, therefore, Lui < if a is sufficiently large. Combining, Lui < in f2 if 
a is sufficiently large. 

It remains to check the boundary condition, 

making use of (H2) and < e < 1. The proof is completed. □ 

Notice that Ui tends to zero uniformly in any compact subset of {x & Q : 
m{x) < ami}. On the other hand, fix any compact subset K of Q\ DJl, 

K <Z [x E Q : m{x) < c^rrii}, 

if we take c < 1 sufficiently close to 1, since all local maximum points of m{x) 
are strict. Therefore, in this case. Theorem 11.51 is a consequence of Lemma 
Oand Lemma O 

Case (6): When m(x) = mi > on 971 and m < at some of its critical 
points. 
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Lemma 2.8. Assume mix) satisfies (H2), (H3) and (H4), and that mix) = 
mi > on SOT. For each c < 1 close to 1, there exists, for all a large, an 
upper solution U2 > in the sense of Definition \2.1\ such that 

,ea{m(x)-crm) ^f^^^ ^(^^^ > q 



U2[x) < 



5 

^a{m{x)~k) when mix) < 0, 



where 0<e<l,A;>0 are appropriately chosen constants independent of a. 

Notice that m. {x E VL : m{x) < ami}, U2 — > as a — )■ oo. We see that in 
this case, Theorem 11.51 follows as before from Lemma [2.81 and Lemma [2.21 

Proof of Lemma \2.8\ . Given c < 1 , let 

9}to = {strict local maximum points xq of m(x) s.t. m(xo) = 0} 
Ai = The union of all connected components of {x G i7 : m{x) > —So} 
not intersecting TIq 

where < 5o < — |max{m(a;o) : x E Q s.t. Vm(xo) = and m{xo) < 0} is 
chosen small enough so that each connected component oi {x E Q : m{x) > 
—So} intersecting TIq lies in {x E Q : m{x) < 0}. This is possible since all 
local maxima are strict. And < e < 1 is chosen to satisfy 

e < (14) 
ami + do 

< k < ecmi. (15) 



k is chosen such that 
Set 



in {x E : m(x) > 0} 
U2= { (j)o in f2 \ Ai 

min{0o5 0i} in Ai \ {x G i7 : m{x) > 0}. 

As before, L0i < in Ai for all a large. On the other hand, by a direct 
computation. 




L(j)o = 0o(^ ~ 0o) < on {x E : m{x) < 0}. 
Hence, Lu2 < for all a large, whenever it is C^. Also, the boundary condi- 

^2 _ ^TTj dm. 



tion ^ — au2^ > is satisfied on dil whenever it is well-defined. 



To see that U2 is an upper solution in the sense of Definition 12.11 it remains 
to show the continuity of U2 and f|T3|) . To this end, it suffices to check the 
following: 
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(i) 01 > 00 in {x eQ : m{x) = -6o}f]d{Ai \{x eQ : m{x) > 0}); 

(ii) 01 < 00 in {x G : m{x) = 0} H d{Ai \{x E^l: m{x) > 0}). 
More precisely, 

(i) : When m{x) = -6o, by (JH]), 

Hence, U2 = 0o in a neighborhood of {x E Q : m{x) = —60} f]d{Ai \ {x E 
n : m{x) > 0}). 

(ii) : When m{x) = 0. bv ffTSj) . 

^ea{m{x)—cmi) ^—eacmi ^ g"'^'^ ga(m(a::)— fc) 

Hence, U2 = 0i in a neighborhood of {x E Q : mix) = 0}f] d{Ai \ {x E Q : 
m{x) > 0}). 

(Notice that 0j are strictly increasing functions of m{x). Hence (possibly 
making 60 smaller) the non-differentiable regions of U2 are regular level sur- 
faces of m{x) by the implicit function theorem.) □ 

Case (c): When m{x) has two distinct values < mi < m2 on 971 and m < 
at some of its critical points. 

We first decompose Q according to the value of m(x). Write Wl = Wli [J 9Jt2, 
where OJtj = {xo G Tl : m(xo) = rrii}, i = 1,2. And define 

dJlo = {strict local maximum points Xq of m{x) s.t. m{xo) = 0}, 



which is possibly empty. Given any c < 1 close to 1, define 



ri = 


{x eQ : m{x) > 0} 














Ai = 


The union of all connected components 
not intersecting Tlo 


of {x 


e n : 


: mi 


[x) 


> 


-5o} 


r2 = 


The union of all connected components 
not intersecting Tli 


of {x 


e n : 


: mi 


[x) 


> 


cnii} 


A2 = 


The union of all connected components 
not intersecting Tli 


of {x 


e n ; 


: m[ 


[x) 


> 


c^rrii} 



where Sq is chosen as in proof of Lemma 12.81 We have a partition: 

= (f] \ Ai) U (Ai \ Ti) U (Fi \ A2) U (A2 \ T2) U 12. 



12 



Lemma 2.9. Given mix) satisfying (H2), (H3) and (H4), and that mix) 
attains exactly two distinct values < mi < 1712 on 971. For each c < 1 close 
to 1, for all a large, there exists an upper solution U3 > in the sense of 
Definition \2.1\ such that 



Usix) < 



^€ia{m{x)—CTni ) 
^e2a(m{x)—cm2) 



in Q\Ai 
in Ai \ A2, 
in A2, 



where < < 1, > are appropriately chosen constants independent of a. 
Notice that in {x G A2 : m(x) < 07712} [J{x G i7 \ A2 : m{x) < ami}, 

M3 — )• as a — )■ 00. 

We see that in the case m{x) having two distinct values mi < m2 on VJl, 
Theorem 11.51 follows as before from Lemma 12.91 and Lemma 12.21 

Proof of LemmalEM Let 0o := e"^"^^^)-'^) and 0i := e^'"(™(^)-'=™') {i = 1,2), 
where < ei < 1 is chosen to satisfy 



5o 



ami + 5o ' 

A; > and < £2 < 1 are chosen such that 

< k < eicmi. 



ei[c-mi - ami 

< 62 < mmj — 

c^nii — cm2 



!}• 



(16) 

(17) 
(18) 



We can now define M3. 



Us 



It can then be proved as before that 
< in 17 and 



00 


in \ Ai 


01 


in Fi \ A2 


< 02 


in 12 


min{0o, 


0i} inAi\ri 


^ min{0i. 


02} in A2 \ r2 



du2 _ dm 

— au2^T— > on aiZ 



whenever they are defined. It remains to show the continuity of U3, as well 
as ffT3ll. It suffices to show: 
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11 

(ii 
[iii 



1 



00 < 01 m{xeQ: m{x) = -Sq} f] d{Ai \ Ti); 

00 > 01 m{xen: m{x) = 0} f] d{Ai \ Ti); 

01 < 02 in {x en-. m{x) = c^rui} f] d{A2 \ T2); 
01 > 02 in {x E : m{x) = cmi} n^(^2 \ r2). 

ii) can be verified following similar lines as in proof of Lemma [221 using 



M and ([T 



(iii) : When m{x) = c^rrii, by (|T8|) 

(iv) : When m{x) = crrii 

^eia{m{x)-cmi) — I ^ ^e2ac(mi-m2) _ ^t2a{m{x)-cm2) q; > 

□ 



Hence, Theorem 11.51 is proved for the cases when m[x) attains 1 or 2 values 
on 971. □ 



The proof of Theorem 11.61 is a modification of the proof in [TT] , overcoming 
the difficulty caused by the local minimum and saddle points of m{x). We 
start with the following lemma. 

Lemma 2.10. With the assumption of Theorem \1.6[ there exists C > such 
that 

u{x) < Ce"("(^'-"^) for all X en and all a large. (19) 
where mi is the unique value of m{x) on 9Jl. 
Proof. Consider w = e^~°'~^'^^™'^^^u{x). Then in f2, w satisfies 

Aw + {a- 2e) Vm ■ Vw - {e{a - e) | Vm|^ + eAm + u- m}w = (20) 
Let z* = z*{a) G be such that w{z*) = max^yw. Then, for x & Q, 

u{x) < ^,(^*)e(-"+^)(™(^*)-'»W). (21) 

We notice that on dQ, 

|^ = e(--M^)(|^ + (-a + e).5^) 
ou ou ou 

ou ou 

du ~ 
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Therefore by the maximum principle, no matter z* G dVt or 'Ww{z*) = 
and Aw{z*) < 0. Hence, by (^U^ 

e{a — e)|Vmp + eAm + u < m at x = z*, (22) 

and 

u{z*) < m{z*) - eAm{z*). (23) 

Now take e = maxj-gj-^^^^}, with the maximum taken over all positive 
saddle points and local minimum points Xq of m{x) such that m{xo) > 0. 
(Take e = 1 if it is an empty set.) Notice that e > by (H4). Then by (122|) . 
we have 

e{a — e)|Vmp < m(z*) — eAm < |m|oo + e|Am|oo, 
which implies that |Vm(2;*)| — )• as a — > oo. Thus, 

dist(z*, {xen-. \Vm{x)\ = 0}) ^ 0. 

Next, we claim that in fact we have dist(z*, OJt) — )■ 0. 

Assume to the contrary that there exists — )■ oo, such that z*{ak) — )■ Xq as 
/c — 7- OO where xq is a saddle point or a minimum point. Then by f l23|) and 
the choice of e, 

< u{z*) < mlz*) — eAm{z*) — i- m(xo) — eAm(xo) < 0, 

which is a contradiction. Therefore, dist{z* ,DJl) — 0. Recalling that m{x) = 
mi on 9Jl, we deduce that there exists C > such that 

mi — m{z*) < C|Vm(2;*)|^, for all a large, 

since the inequality holds in a neighborhood of Wl, where z* eventually enters. 
Hence by (l22l) again, 

(a - e)(mi - m(2*)) < C{a - e)\Vm{z*)\^ < C(^^^ - Am{z*)). 
Therefore, 

(a - e)(mi - m(2*)) < C(— + IIAmlL) (24) 

e 

And for every x E Q, from fl?I]) . 

= ■u(2*)e'^'^™i~'"*^^)''"'''-"~'^^*^'"i~™*^^*''^ 

< (mi + e||Am||oo)e''l'"l°°+^(^+"^""°°\ 
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by fl2^ and (^^. Since the right hand side is a constant independent of x 
and a, ffTI?]) is proved. 

□ 

Proof of Theorem \l.(A From (|T9l) . we see that for all p > 1, m — in as 
a — )■ oo. For each xq G fix a neighborhood il(xo) of Xq, by (IT^ . 

where Ci = ||D'^m||oo/6. Denote M(xo,a) = supj^(xo)M, which is attained 



in Bjii^[xq) for i? sufficiently large, and all large a (by Theorem 11.41 and 
Lemma [2.10p . Define 

u{xq + 

Then sup Wa = 1 in •\/a(il(xo) — xq), and 

for all a large and in {y G Mat : Xq + y/y/a G fi, |?/| < -^^^^}, where 
Ai < ■ ■ ■ < Aat < are the eigenvalues of L'^m(xo). 

To prove (jl]), by Lemma [2.101 and the fact that M(xo,a) is bounded, it 
suffices to show that for each xq G £0T 



Wa{y) qIv"^ D^rn{xo)y -^^ gvcry compact subset of M^, and 
M(xo,a) ^ 2^/2^(xo), 

as a oo. Wa satisfies /^yWa + ■ V^PVa + QW^o = 0, where 

^ = ^(a, y) = -i/a ■ WxmixQ + ^) , 

A/a' 

and 

y u{xq + -^) — m(xo + 



(25) 



'a a 
The boundedness of u (by (fT9|) ) implies that 

lim '^{a,y) = -y^D'^m{xo), lim Q{a,y) = -A^m(xo), 

a— ^oo a—^oo 

uniformly in any compact subset of M^. Hence by elliptic estimates (see 
[16]), using the fact that for each compact subset K in M^, Wa is bounded 
in W^K) for p G (1, oo] and all large a, after passing to a subsequence if 
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necessary, as a — ?■ oo, Wa converges to some function W* uniformly in any 
compact subset of M^, and W* must satisfy 

AyW* - yD^m{xo)VyW* - Am{xo)W* =0 in M^, 
supj^^ W*{y) = l, < W*{y) < Cei^^^'"^^^")?' \/y G M^. ^ ' 

Now we invoke the following lemma, the proof of which makes use of a 
Liouville-type result due to [17] which is formulated differently in [15], and 
will be included in the Appendix C for completeness. 

Lemma 2.11. If W* G Wl^^{R^) satisfies then W* = ei2'^^''"(^o)2^. 

The uniqueness of the limit implies that 



lim Wa{y) = ^ -^My uniformly in any compact subset of . (27) 



That W* attains its strict maximum at the origin and f[T9l) implies that 

lim ^r^j^ = W*iO) = 1. (28) 
To show the second part of (!25|) . it remains to calculate lim u{xo). In [TT] 

a— >oo 

it was accomplished when m as a single peak via a "global" argument. Here 
we devise a "local" argument near each xq G DJl. 

Lemma 2.12. For each xq G Wl, liminf ^(xo) > 2^^^mi. 



Proof. By following the proof of Theorem 11.41 with the same choice of test 
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function ■?/' and open sets f/j, we have for each r] > 0, 
< hm inf 



< hm inf 

a— >oo 



■BH/v^(a;o) 



; dx Jj 



/c/2 e^t"^-™!] (ix 



+ 



1^ dx 



J^^ e»[^-"^i] dx 



hm 



< hm inf 

a—^oo 



+ 



< hm inf 

a— >oo 



- "^(xo) 

a[m{xo + ^)-mi] + ^y'^ D^m{xo)y 



dy 



(0) 



-cils/P dy' 



dy 



m(xo) 



iV 

<(l+r7) hminf'u(xo) (2~^ + r^) + 77 — m(xo) 

The third inequahty fohows from f l27|) . f l28|l and the Lebesgue Dominated 
Convergence. In the fourth inequahty, we apphed the change of coordinates 
X = xo + ^ and that there exists ci, C2 > such that ci|?/p < mi — m{x) < 
(which are consequences of the nondegeneracy of m). The last hne 
foUows by taking R > sufficiently large and that 

y 1 

lim a[m{xo H 1=) — mi] = -y'^ D'^m{xo)y 

a-i-oo yja 2 

uniformly in compact subsets of R^. Finally, the lemma is proved by letting 
r]^0+ □ 

Next, we claim that 



Claim 2.13. lim V / e'^y^^'"'^'"'^^ dy \u{xof - 2^^^ 



miulxo 
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Proof of Claim l2JR Integrate ( ITUj) over Q, we have 



= {u — um) dx 
'n 



+ j + j ^ {u^ — um) dx 

UiDi-Bfl/v^(a;o) -'Un)t-Bro(a^o)\-Bfl/y^(a:o) J n\UsnBro(xo) 



E 

+ 0(e- 



/ (m^ — um) dx + C I 

JBa/^{xo) Jb 



■■ya\ 



by Theorem II .51 and Lemma l2.10[ Muhiply by a 2 and changing coordinates 
X = Xq + we see that 

0= V / {u^-um){xQ + ^)dy + 0{l e-'^^^^" dy)+0{a^ e-'"'). 

By (1271) and (1251) . for each > large, there exists ao such that for any 

a > ao, 



dy + o{l) 



+ 0{ 



dy) 



m\BR{o) 



dy + 0(1) 



e-^'^y^'dy). 



where hm o(l) = 0. Now take a — )■ 00 and then R — )■ 00, we have the 



desired result. 



□ 



Lemma 12.121 and Claim 12.131 implies the second part of fl25l) . This concludes 

□ 



the proof of Theorem 11.61 



3 Proof of Theorem 11.9 



As before, assume for simphcity di = 1. 
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Proof of Theorem \1.9\ . Notice that (f/^, Va) satisfies 

V ■ {\/U - aUVm) + U{m - U) = UV > in VL, 

d2AV + Vim -V) =UV >0 in (29) 
- aU^ = |K = on dn. 

By method of upper and lower solutions, < f/^ < m and Q < Va < Od2- 
([7]) follows from the same argument as in proof of Theorem ll.4[ using the 
inequality < ^^2- That Ua converges to away from the positive local 
maximum points of m{x) follows from the corresponding property of u. 

Now, assume m = mi on the set of its local maximum points. 

Lemma 3.1. // (H2), (H3) and (H4) hold, and m{x) is constant on its 
local maximum points, then there exists C2 > such that 

Ua{x) < Cse"^"'^^)-'^^) for all x e Q and all a large. 

Lemma [3.11 follows from Lemma [2.101 and the fact that < [/q, < u. 

For some large, JqItti — 6*2 e""'-'"^^^"™^^] > and by a claim on P. 498 in 
[H], there exists a positive solution Vq of 

d2AVo + Vo{m - (:72e"«(™(^)-™i) - V"o) = in Q, 
§a = on dQ. 

av 

then for all a > ao, 

AVq + 1/o(m - f/a - K)) > in Q, 

^ =0 on an. 

Therefore, Vq is a lower solution of the second equation of f l29|) for Va, and, 
6d2 >Vo,>Vo>0 for all a > a^. (30) 



By Lemma 13. 1[ ?7q, — )■ in for any p > 1. By second equation in 
(150]) . and elliptic estimates and uniqueness, V 6*^3 weakly in iy^'P(f2) in 
any p > 1 hence strongly in C^'^{VL) for any (3 G (0, 1). This proves (j8]). 

Fix xo e 9Jt and let W^{y) = ^°i7(+y , where M(xo,a) = sup^„^^(,^,) 
for some small tq > 0. (M(xo,a) is independent of the choice of tq by ([7]) 
and Lemma [3.1[ ) As in proof of Theorem 11.61 notice that Wci{y) — > W*{y) 
as a — 00 uniformly for y in compact sets in where W* satisfies 



\W* - yD'^m{xo)VyW* - Am(xo)IV* =0 in 



nN 
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Also similar as in proof of Theorem 11.6^ 

lim WM = W*{y) = ei?'^^'™(^'o)?^ on compact sets in (31) 

a— >-oo 



and lim — = 1. Now, by arguments in the proof of Theorem II ■4[ we 

a^oo M[Xi, a) 

have 



m] 



liminf ^ — > 0. 



a^oo / gam 

J U ~ 



1U2 

Then Lemma 13.11 ([H]) and (ISTj) implies, for each xq G 971, 

liminf f/«(a;o) > 2^/\m, - e,,{xo)) (32) 

a— >oo 

By integrating the first equation of (l29|l over Q, we have J^Ua{rn — Ua — 
Va)dx = 0. And by similar arguments in proving Claim [TT^ we have 

0= lim V / e5^"^M^»)^rfi/[t/„(xo)^-2^/2(mi-^^,,(xo))f/„(xo)]. (33) 
Finally, lim U^{xo) = 2^/\mi - ^^^(xo)) follows from ([32D and □ 



4 Concluding Remarks 

In this paper, the existence of concentration phenomena in the globally sta- 
ble steady state u{x) of ([2]) is proved for m{x) which has finitely many local 
maximum points. Furthermore, the concentration set is shown to be the set 
of positive local maximum points of m{x). The situation when m(a;) con- 
tains local maximums that are not strict is however, completely open. It is 
possible that u would concentrate on some higher dimensional sets. 

In this paper, the limiting profile is obtained in the special case when the 
resource function m has equal peaks. Based on the estimates established in 
this paper, a special method is introduced to determine the limiting profile 
for m with peaks of different heights in [18]. However, the method only works 
for = 1. For N > 2, very recently the limiting profile has been found by 
the author. This will be published in a forthcoming paper. 

We learnt recently that in [19], a lower solution for ffTOj) can be constructed at 
each xq E DJt which gives an alternative proof for the existence of peaks on DJt. 
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We also remark that the assumptions on m{x) in {x G : m(a;) < 0} can be 
weakened substantially. In fact, instead of (H2), (H3) and (H4), we only 
need to assume that there exists S > 0, such that the foUowings hold. 

(H2') ^<Oon{xedn: m{x) > -6}. 

(H3') m(a;) has finitely many local maximum points in {x E Q : m{x) > 
—6}, all being strict local maxima and are located in the interior of Q. 

(H4') If xo G satisfies m(xo) > —S and is a local minimum or a saddle 
point of m{x), then Am{xo) > 

Finally, notice that although we have set the diffusion coefficient d,di = 1 
for simplicity, the results proved in this paper hold true for any d, di > 0, as 
stated in Section 1. 

5 Appendix A 

Denote 9d to be the unique positive solution to 

r dAO + e{m-e) =0 in n, 
\ ^ =0 on on. 

The existence part is standard. (See, e.g. P. 498 in [9].) Also, it is known 
that (Prop. 3.16 of [I]) 

lim 9d = m'^, uniformly in f2, (Al) 

d-s>0+ 

where m^{x) := max{m(x), 0}. 

Here we shall prove that if xq € f2 is a positive strict local maximum point of 
m and Am(xo) < 0, then m(xo) — Od{xo) > for all > sufficiently small. 

Remark 5.1 ([20J). When d is not small, there are counter examples showing 
that the conclusion is not true in general for xq G SOT other than the global 
maximum point(s). 

First we show that m(xo) > ^^(^^o)- Assume now to the contrary that for 
some positive strict positive local maximum point xq of m(x), for some se- 
quence c?i — >■ 0, 

9d,{xo) > m(xo) > 0. 

Now, Xq E {x E Q : 9di{x) > m(xo)} for all i. Denote by Ui the connected 
component of {x G f2 : 9di{x) > m(xo)} that contains xq, then Ui ^ ^ and 
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^dfii = Gdi{(^di — rn) > in Ui. i.e. is subharmonic in f/j. Now for di 
sufficiently small, by (Al), Ui is compactly contained in a neighborhood of 
xq. In particular, 6'^- > m{xo) in Ui and Odi{x) = m{xo) on dU. This con- 
tradicts the property of subharmonic functions. Therefore, m{xo) > 9d{xo) 
for all d > sufficiently small. 

Now assume there exists a sequence rf, — )■ such that 6'rf-(xo) = m{xo). We 
claim that 

Claim 5.2. VOdX^o) = for all i sufficiently large. 

Otherwise there exists Xi — )■ xq such that OdX^i) > rn{xo) and a contradiction 
can be reached by previous arguments by choosing a horizontal hyperplane. 

Now since OdX^o) = m{xo), VO^Xxq) = Vm{xo) and A6d, = > Vm(xo), 
there exists Xi — )■ Xq such that Od-{xi) > m{xi). (Since otherwise the mean 
curvature of the surface defined by 9d.i in ]R^+^ at xq, which is a multiple of 
A0d2{xo), would not be not equal to 0.) Now fix a neighborhood Uq of Xq, 
and a (slightly tilted) hyperplane Sj : L(R^,]R) such that 

9dXxi) > T.i{xi) and T.i{x) > m{x) in Uq. 

By (Al), 9di — > m uniformly on dU^ while mina;7o{Sj(x) — m(x)} > c > 
for some constant c independent of i. This implies that there is some Ui 
such that 

r A^rf, =erf,(^^rf^-m)>Oinf/, 
\ 9d.^ > Sj in f/j, Od.^ = Sj on c^t/j 

which again contradicts the fact that Od^ is subharmonic in Ui. 

6 Appendix B 

Here we discuss the proof of the general case of Theorem II. 5[ Recall 

DXt = { positive strict local maximum points of m{x) in Q }. 

By (H3), m{x) has finitely many local maximum points. Let < mi < 
m2 < ■ ■ ■ < be the distinct values of m{x) on 9Jt. Decompose 

no 

1=1 

where Wli = {xq G Tl : m(xo) = nii}. And let 

VJto := {local maximum points xq of m{x) s.t. m{xo) = 0}, 
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which is possibly empty. For each c < 1, close to 1. Define 60 as in the proof 
of Lemma [2.81 Decompose fl according to the value of m(a;): 

Ti ={x G : m{x) > 0} 

Ai =Union of connected components oi {x & Q : m(x) > — 5o} 

not intersecting 9Jlo, 
Fj =Union of connected components of {x E fl : m{x) > cnii^i} 

not intersecting VJti-i, 
Ai =Union of connected components of G : m{x) > c^nii-i} 

not intersecting 

for i = 2, rio- Notice that Aj D Fj D Aj+i 3 Fj+i. Define 

for i = 1, no — 1 
= < e"("'(^)-'=) in n \ Ai 

jj^j^|geia(ni(a;)-cmi)^ gei+ia(ni(x)-cmi+i ) j ^.^^ y "p^^^ 

for i = 1, 77-0 — 1 

where 0<ej< l,/i;>0 are constants chosen such that 

ei < —, < k < eicnii, and 

crrii + do 

< ej+i < mmj^r , 1|, for z = 1, ■ ■ -, Uq - 1. 

Then, we have 

Lemma 6.1. Given m(x) satisfying (H2), (H3) and (H4). For every 
c < 1 sufficiently close to 1 , u > is an upper solution to ( ITOi) according to 
Definition \2.1\ 

The proof of Lemma 16.11 is similar to that of Lemma 12.91 and Lemma 12.81 and 
is omitted. 

Notice that the full statement of Theorem 11.51 follows from the above lemma 
and Lemma [2. 2[ 
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7 Appendix C 



Next, we shall prove Lemma 12.111 We first state and prove the following 
Liouville-type theorem which is due to |17j, following the formulation in |15] . 



Theorem 7.1. Let a G L^j,(M^) be a positive function. Assume that $ G 
(M^) satisfies in the weak sense 

$ div{a^V<^) >0 in , 

and for some C > and every R > 1, 



Jbr{0) 



Then ^ is a constant. 



(CI) 
(C2) 



Proof of Theorem \7. 1\ From (1C1|) we deduce, for any smooth function tp, 

div($^V2V$) > ^V^l V$|^ + 2^ipa^V'4; ■ V^. (C3) 

Let C be a C°° function on [0, oo) with < C{t) < 1 and C(t) = 1 for 
< t < 1, Cit) = for t > 2. For > and s G set Cnix) = Ci\x\/R). 
Taking ip = in (\C3\i and integrating over M^, we find, by the divergence 
theorem, 



CV|V$prfx < 2 



< 2 



ni/2r 



a^^'^lVCnl'dx 



'R<\x\<2R 

By (IC2p and the definition of (ji, we can find Ci > such that 



1/2 



Therefore 

[ C],a^\V^\^dx<2^, \ [ CW\V^\^dx 

Jr'^ IJr<\x\<2R 

This implies that 



1/2 



(C4) 



Cla^lV^l^dx < 4Ci, 



and hence, letting — )■ cxd in (IC4p we obtain 



a^WWdx = 0. 



This implies |V$| = a.e. Hence $ is a constant. 



□ 
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Proof of Lemma \2.11[ Given W* satisfying (12^ . we want to show that W* = 

First we make the transformation W* = e~3y^^^™'(^o)2'$. By (12^ . we see 
that $ satisfies 

div(e52^^^'™(^o)?'V$) =0 in M^, 
< $ < i^'3e-^2'^^''^(^o)2', 

SUpKiV $(|/)e52^^^''™(^«)?' = 1. 

It remains to show that $ is a constant. By Theorem 17. it suffices to show 
that for some C > and every R> 1, 

f e^y^^^'^^'^'^y^^dx < CR". (C5) 

Jbr{0) 

By noticing that the integrand can be dominated by 

we have immediately that flC5l) is true. Hence the theorem is proved. □ 
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